Magnetic frustration on the diamond lattice of the A-site magnetic spinels 
CoAl2-aGa^.04: The role of lattice expansion and site disorder 
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The spinels C0B2O4 with magnetic Co^"^ ions on the diamond lattice A site can be frustrated 
because of competing near-neighbor (Ji) and next-near neighbor (J2) interactions. Here we describe 
attempts to understand these interactions by substitution on the non-magnetic B-site. The system 
we employ is CoAl2-i,Ga3;04, where Al is systematically replaced by the larger Ga, ostensibly on 
the B site. Ga substitution has the effect of expanding the lattice and pushes Co atoms on the 
A-site further away from one another weakening the magnetic interactions. We also find, however, 
that Ga distributes between the B and the A site in a concentration dependent manner displacing 
an increasing amount of Co from the A site for larger values of x. This site mixing, confirmed 
by powder neutron diffraction studies carried out at room temperature, affects magnetic proper- 
ties very significantly and changes the nature of the ground state. We discuss the role that both 
structural changes play in changing the degree of magnetic frustration on the diamond lattice. We 
also use classical Monte-Carlo modeling of the magnetic coupling to illustrate the complexity of the 
interactions that arises from site mixing. 



PACS numbers: 75.50.Ee,75.40.Mg, 

I. INTRODUCTION 

The presence of magnetic frustration in the solid state 
is known to give rise to multiple, nearly de gener ate and 
often non-coUinear magnetic ground states .EEEl Interest 
in the competition between these ground states has at- 
tracted considerable attention in recent years with the 
resurgent field of multiferroics: systems with spiral spin 
ordering are k nown to couple spin and lattice degrees 
of freedorrPI^ and in specific cases, to even give rise to 
spontaneous polarization.^^ 

One highly studied class of materials that has been 
shown to exhibit this type of non-coUinear ordering re- 
sulting from magnetic frustration is the spinel family.l^Eni 
Spinels have the general formula AB2X4, and in so-called 
normal spinels, the A site is divalent and tetrahedrally 
coordinated by the anion X, while the B site is trivalent 
and octahedrally coordinated by the anion X (typically 
O or S). Both sites can accommodate magnetic cations 
allowing for a wide range of magnetic properties. It is 
important to note that these sites are often not well or- 
dered and as a consequence it is frequently appropriate 
to write the formula as {Ai^s^s){.^s'B2-s)^4 where 6 is 
referred to as the inversion parameter and can be any 



value <6 <1. 

The spinel structure can be viewed as two interpen- 
etrating sublattices with B atoms forming a pyrochlore 
lattice while the A atoms constitute a diamond lattice.!^ 
It is well known that there are multiple magnetic inter- 
actions in these systems. Here we focus mostly on the 
nearest neighbor (Ji) and next nearest neighbor (J2) 
exchange interactions on the diamond lattice using the 
compound C0AI2O4 as a starting point (Fig.[T|). Since 
the diamond lattice can be broken down further into two 
interpenetrating face-centered cubic (fee) sublattices, we 
can envisage Ji as coupling the two fee sublattices to 
each other and J2 coupling nearest neighbors within each 
sublattice-li^l 

In a recent work Tristan et al. demonstrated that 
normal spinels, 5 <C 1, with magnetic cations on the 
A site and non-magnetic cations on the B site can ex- 
hibit strong magnetic frustration .^^^ These A-site mag- 
netic spinels have been studied in the past by RotH^ 
and Blasse.ES Recent modelling studies by Bergman et 
al. suggested that the magnitude of the ratio J2/J1 can 
strongly influence the dominant magnetic ground state 
in these systems-^^In the limit that J2/J1 ~ 0, the mag- 
netic ordering is a simple Neel state (magnetically speak- 
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FIG. 1: (Color online) Illustration of the superexchange path- 
way for magnetic interactions between neighboring A-site 
atoms in the AB2O4 spinel structure. A atoms are blue (dark 
grey), B are light grey (light greay) and O are the orange 
(grey) spheres. Note that the near-neighbor Ji coupling pro- 
ceeds from each A atom to one of it's 4 near-neighbor A 
atoms through this 6-fold degenerate pathway. The next- 
near-neighbor J2 coupling to one the 12 next-near neighbors 
is doubly degenerate and proceeds through an identical path- 
way. 

ing, a diamond to zinc blende transition). As the ratio 
becomes larger, the nature of the magnetic ground state 
can change from Neel to a complex spiral pattern with 
the spirals oriented parallel to (111) planes. Tristan et 
oZ.I^ have also recently examined solid solutions between 
C0AI2O4 and C03O4 to investigate the effect of the B 
site cation in the superexchange pathways that connects 
neighboring A sites. 

These recent studies have encouraged us to attempt to 
alter this ratio of J2/J1 in these A-site magnetic spinels 
by controlling interatomic spacings through appropriate 
substitution on the B site of the spinel structure. We 
use Ga substitution ostensibly on the B site of spinels 
CoAl2_a;Gaa;04 as a means of separating the Co from 
one another, given that Ga''"'" is significantly larger than 
AP"*" (the respective 6-coordinate radii are 0.62 A and 
0.535 A) We find that while the effect of such substi- 
tution is indeed to steadily increase the distance between 
neighboring Co atoms, there is the additional effect that 
site inversion in the structure steadily increases simul- 
taneously. While separating the effects of these simul- 
taneous changes in the structure is dificult, we discuss 
the role each change plays in the context of the altered 
frustration parameter. 

II. EXPERIMENTAL DETAILS 

Polycrystalline samples of the compounds were 
prepared using ceramic routes. Cobalt oxalate 
(CoC204-2H20) was mixed with stoichiometric amounts 



of Ga203 and AI2O3 and intimately ground with ethanol 
in an agate mortar. The powders were then pressed into 
13 mm pellets and fired in air in alumina crucibles at 
800°C for 24 h. The pellets were then reground, pressed 
again into pellets, and fired at 1200°C (1000°C for x=0.0) 
for 12 h. In order to obtain equilibrated samples, all pel- 
lets were annealed by heating to 700°C for 12 h, cooling 
to 400°C at a rate of 3°Cmin-i, soaking at 400°C for 
120 \^ followed by cooling in the furnace to the room 
temperature. For all heat treatments, pellets were placed 
on a bed of powder with the same stoichiometry to min- 
imize reaction with the crucible. 

X-ray diffraction patterns were obtained using CuKa 
radiation on a Philips XPERT MPD diffractometer op- 
erated at 45 kV and 40 mA. Phase purity was determined 
by refining the patterns using the Rietveld method as im- 
plemented in the xnd Rietveld code.'^ Neutron diffrac- 
tion data were collected on the neutron powder diffrac- 
tometer (NPDF) at the Lujan Center at Los Alamos 
National Laboratory at room temperature, on samples 
sealed in vanadium cans .'^Sl Neutron diffraction data were 
refined using the Rietveld method as implemented in the 
EXPGUI-GSAS software suite.l^iJSll Local structures as ob- 
tained from pair distribution function (PDF) analysis of 
the total neutron scattering were extracted from the total 
scattering data using the program PDFgetl^P^ and with 
a maximum momentum transfer Qmax = 40 A~^. The 
obtained PDFs were analyzed using the PDFgui soft- 
ware package.I^ DC magnetization was measured using 
a Quantum Design MPMS 5XL SQUID magnetometer. 

III. COMPUTATIONAL DETAILS 

Classical Monte Carlo simulations of the magnetic be- 
havior of the system were performed using the ALPS 
project's spinmc application.!^ A custom lattice with pe- 
riodic boundary conditions was generated for each sim- 
ulation run in order to allow disorder averaging. The 
lattice generation code first randomly selected an A site 
and then moved its spin to a B site, also chosen at random 
and independent of the position of the A site. The simu- 
lation was then run using a Heisenberg Hamiltonian with 
J I and J2 bonds as described above, and with Ji impu- 
rity bonds connecting spins on B sites with their nearest- 
neighbor occupied A sites. Finally, because the typical 
error bars obtained from each Monte Carlo susceptibility 
simulation were negligible (about 1 x 10""* relative er- 
ror), the error bars plotted represent only variation due 
to the presence of disorder as obtained by averaging over 
several runs. 

A numerically calculated inverse susceptibility curve 
was fit to the experimental curves as follows: first, using 
the fact that the Curie- Weiss temperature for the J1-J2- 
Ji model with an inversion parameter, 5, is given by 
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FIG. 2: Highest d-spacing bank of powder neutron diffrac- 
tion data from the different samples of CoAl2_iGaa;04 col- 
lected at room tempearture. Circles are data and solid lines 
are Rietveld fits to the spinel phase [space group Fd3m, A 
(|, |, |), B (|, |, |), O {u,u,u) ]. Asterisks denote locations 
where contributions to the scattering from the AI2O3 impu- 
rity arise. 



one can calculate the value of Ji necessary to match the 
experimentally measured Curie- Weiss temperature. This 
value of Ji is then used to scale the temperature by a 
factor of Ji S{S + 1) and the inverse susceptibility by 
a factor of {ks Ji S{S + I)) / {^l^^'^gNA)- Such a fitting 
procedure always guarantees that the high temperature 
behavior of the simulation data is in exact agreement 
with that of the experimental data. 



IV. RESULTS 

A. Structure 

Rietveld analysis of the highest c?-spacing bank (from 
four banks of data for each sample) of time-of-flight neu- 
tron diffraction data are displayed in Fig.|2]for the differ- 
ent spinel samples. There were no peaks in the diffraction 
data that could not be assigned to either the spinel phase 
or a small AI2O3 impurity in the x = 0.0 sample, indi- 
cating no magnetic impurities and that a complete solid 
solution is achieved for all values of x. The oxygen stoi- 
chiometry refined within error to the correct stoichiomet- 
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FIG. 3; Evolution of different refined structural parameters 
with X for the different CoAl2_a;Ga2;04 samples obtained from 
time-of-flight neutron diffraction, (a) The cell parameter, 
showing a nearly increase in size. The line connects end 
members to illustrate the Vegard law. (b) Occupancies of 
the different cations on the A site, (c) Internal parameter 
u indicative of the oxygen position, u — ^ means the BOg 
octahedra are perfectly regular. 



ric value excluding the possibility that any of the Co^+ 
was oxidized to Co'^"'". Given the close-packed nature of 
the oxygen lattice this was not surprising. From the Ri- 
etveld analysis, evolution of the lattice parameter, site 
occupancy of the cations on the A site, and the inter- 
nal structural parameter (the u of oxygen) are presented 
in the different panels of Fig. [3] It is seen that substi- 
tution of the larger Ga'^"'" for Al^"*" results in the spinel 
unit cell edge increasing from 8.1 A to 8.3 A. The Vegard 
law is not strictly followed, and for all intermediate x val- 
ues, the cell parameter is slightly reduced from the values 
suggested by a weighted average of the end-members. 

Panel (b) of Fig.js] shows the results of allowing all 
ions, Al, Ga, and Co to distribute themselves between 
the A and B sites in the refinements, with the con- 
straint that the total amounts of the different atoms were 
as dictated by the starting stoichiometry. The process 
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of achieving this refinement involved declaring on each 
of the A and B sites, two separate Co atoms, one of 
which exchanged with Al and the other with Ga. A re- 
lated procedure has been described by Joubert et alW^ 
Due to the nature of multiple constraints in the refine- 
ments, errors on the occupancies are likely to be un- 
derestimated. It is noted from Fig.jsj^b) that as larger 
amounts of Ga are substituted into the system an in- 
creasing displacement of Co from the tetrahedral site to 
the octahedral site occurs. While it might be expected 
that the Ga'^"'" ions would prefer the octahedral site based 
upon its large size, the amount of Ga^"*" found on the 
smaller tetrahedral site steadily increases across the sub- 
stitution series. In contrast Al'^+ ions are found to re- 
main mostly normal with a tetrahedral site occupancy 
never exceeding 20%. The site preference of Ga^+ for 
the tetrahedral site is in agreement with the tendency 
of Ga^^ to adopt sp'^ hybridization. In fact, from 
electrostatic arguments Millei^^ has determined that the 
relative octahedral site preference of Al^+ and Ga'^"'" are 
respectively — 2.5kcalmol~^ and — 15.4kcalmol~^ with 
the larger number indicative of the greater octahedral 
preference. In contrast, Co^+ has a site preference en- 
ergy of — 110.5 kcalmol"^. Nakatsuka et al^have con- 
ducted a recent study on the energetics of different local 
bonding configurations and found that replacing the rel- 
atively large Co^+ {rtet — 0.58 A) on the tetrahedral site 
with Al^^^ {rtet = 0.39 A) results in abnormally long bond 
lengths which is not favored. This effect is not as pro- 
nounced when Ga'^''" moves to the tetrahedral site given 
the larger radius {rtet = 0.47 A) and that the end mem- 
ber CoGa204 is commonly found to be almost completely 
inverted.'^ Site mixing also serves to explain the devia- 
tion from the Vegard law that is observed in Fig.jsja). 

In Fig.jsjc), the internal parameter u reflecting the po- 
sition of oxygen is displayed for the different compounds 
in the series. It is seen that with increasing x, this value 
progressively decreases. Hill et al.^'^ have pointed out 
that in normal spinels (without any inversion) this pa- 
rameter u depends on the ratio R of the octahedral to 
tetrahedral bond lengths according to: 

i?V4 - 2/3 + (lli?V48 - 1/18)1/2 
2W^2 

The expected values for ordered C0AI2O4 and CoGa204 
using the appropriate ionic radiP^ would be respectively 
0.265 and 0.261. The observed trend of a decreasing u 
agrees with this expectation for the normal end-members. 
However, we see from Fig.jsjc) that while C0AI2O4, with 
very small inversion has an experimentally determined u 
value very close to what is calculated from equation[2] 
the value determined for CoGa204 {u = 0.258) is sig- 
nificantly smaller than suggested by equation[2j Again, 
we believe this discrepancy arises because of the growing 
inversion in the compounds as x increases. 

The pair distribution functions extracted from total 
scattering neutron diffraction are shown for each compo- 
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FIG. 4: (a) Neutron pair distribution functions (PDF) col- 
lected at room temperature for the different CoAl2-a;Gai:04 
samples with x = 0.0 at the top of the panel and x = 2.0 
at the bottom. Experimental data are filled grey circles and 
the lines are fits obtained by refining data against the aver- 
age spinel structure, (b) Experimental PDFs of the x = 0.4 
sample (filled grey circles) compared with the refinement us- 
ing the average structure (labeled "fit") and compared with 
a simulation (labeled "simulation") prepared by taking the 
weighted average of experimental PDFs of the x — 0.0 sam- 
ple (80% weighting) and x = 2.0 (20% weighting). 

sition in Fig.|4|^a) along with fits (shown as lines) obtained 
by refining the average spinel unit cells with the appropri- 
ate site inversion. Several characteristics can be observed 
directly. The unit cell expansion is apparent in the move- 
ment of atom-atom peaks to longer r with increasing x. 
The first atom-atom peak become increasingly broad as 
X increases, perhaps as a result of increased site mixing, 
in turn resulting in a distribution of cation-oxygen dis- 
tances. Furthermore, features in the intermediate com- 
positions can be linked to the C0AI2O4 and CoGa204 
end- member PDFs, with an increasing resemblance to 
the CoGa204 PDF as x increases. 

Refinements in real space were carried out starting 
with the average structure results. A stable refinement 
of the site occupancies could not be achieved and val- 
ues were therefore fixed to the results obtained from the 
Rietveld analysis of the neutron scattering data. The pa- 
rameters refined in the PDF analysis include lattice pa- 
rameters, isotropic atomic displacement parameters for 
each atomic species, a scale factor, and quadratic peak 
sharpening. The refinements yielded < 1 ^^'^ Rwp 
less than 15% for r-ranges from 1 to 20 A. From the fits 
of the PDF displayed in Fig.|4|^a), it is seen that of all 
the samples, only in the data corresponding to the x = 
0.0 sample are all the peaks well described by the model. 
This is keeping with the fact that only in the x = 0.0 
sample docs the average spinel structure accurately de- 
scribe all the local distances, since this structure has the 
lowest inversion and no substitution. 



5 



In Fig.|4]ja), this is shown for a particular sample with 
X = 0.4. The experimental PDF is only very poorly fit in 
the short r region by the average structure model (labeled 
"fit"). Upon close inspection of the a: = 0.4 refinement 
around 3 A and 5 A, which correspond respectively to the 
first and second nearest Co neighbors, there is significant 
splitting of the experimental peaks which is not well rep- 
resented by the average structure. However, a much bet- 
ter description of the a; = 0.4 sample is a stoichiometri- 
cally weighted (80:20) average of the experimental PDFs 
of the end members (x — 0.0 and x — 2.0) rather than 
as a single phase as illustrated in Fig.|4](b). In this, a 
resemblance to systems such as Ini_3,Gaa;As^ is noted 
wherein the alloy compositions locally follow the bonding 
rules of the end-member structures, and alloy PDFs can 
be described using weighted averages of the end mem- 
bers. This result of averaging end- members only holds 
true for distances within a unit cell (approximately 8 A). 
Outside of this range the superposition model begins to 
fail and the average structure proves to be a good model 
of the data. Thus from the local structure analysis we 
can see that our samples have a homogeneous distribu- 
tion of substituted cations and there does not seem to be 
any evidence for local clustering of cations. 



B. Magnetism 

Figure|5] displays the temperature dependence of the 
DC magnetic susceptibility x of the different spinel sam- 
ples on a single scaled plot. For all samples, data be- 
tween 350 K and 400 K were fit with high reliability by 
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the Curie- Weiss formula: x = C'/(r— Qcw), where C is 
the Curie constant, and Qcw is the Curie- Weiss ordering 
temperature. We can recast this formula as: 
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FIG. 5: (Color online) Temperature dependence of the scaled 
zero-field cooled magnetic susceptibility of the different spinel 
samples. Data were acquired under a 100 Oe field. The nature 
of the scaling is described in the text. 
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Evidence of Curie- Weiss behavior at high temperatures is 
seen from the right, higher temperature, panel of Fig. [5] 
For each sample, and at all temperatures above the in- 
dividual Qcw (indicated in Table|l] and approximately 
ranging between -100 K and -40 K) the scaled inverse sus- 
ceptibility (left-hand side of equation|3| is precisely equal 
to T/\Qcw\- AH the samples order at temperatures well 
below |6(7w| as seen in the left-hand low-temperature 
panel of Fig.[5j A curious point to note is that while all 
the samples display a downturn in plots of x W'S- T, it 
is only the x — 0.0 sample, C0AI2O4, that stays largely 
superior to the dotted Curie- Weiss line. The behavior 
oi X — 0.0 is distinctly different from all other samples 
which progressively deviate from the dotted Curie- Weiss 
line at higher temperatures as x increases. The utility 
of plotting the inverse susceptibility of a solid solution in 
the manner shown in Fig. [5] becomes evident in the way 
compounds with uncompensated spins due to inversion 
{x > 0.4) are separated from the antiferromagnetic x = 
0.0 end member^^. A plot of 1/x vs. T does not re- 
veal this since the magnititude of the susceptilities varies 
through the series. 

Table|l] also shows values of fi^f / obtained from the 
Curie constant C for the different spinel samples. The 
expected spin-only value of the magnetic moment for 
tetrahedral Co'^^ is 3.88 hb whereas a value of 5.20 ^ij^ 
is expected for systems with completely unquenched or- 
bital contribution. The values obtained here run between 
4.59 and 4.85 and are therefore sandwiched by 
the limits of the completely quenched and unquenched 
orbital contributions. Measured values from the liter- 
ature are in the range of AA^b to 4.8 /is !^ Cossee 
and van ArkeP^ have argued that for tetrahedral Co^+ , 
the proximity of a low-lying excited spin state adds a 
temperature-independent term to the Curie- Weiss law 
and that after making such a correction, magnetic mo- 
ment is close to 4.4 ^b- We expect C0AI2O4, with a very 
small degree of inversion, to display a Me// value which 
is close to AAfiB- We find that for x ^ in the sub- 
stitution series, the value of Heff is always larger than 
for X — 0.0, which we attribute to the increasing amount 
of octahedral Co^+ which has a larger orbital contribu- 
tion to the effective moment and correspondingly has an 
experimental range of fieff from 4.7/iB to 5.2 fXB^ 

Figure[6] shows in closer detail, the temperature depen- 
dence of the magnetic susceptibilities of the title spinel 
compounds at low temperatures. All the compounds dis- 
play splitting of the FC and ZFC data between 4 K and 
12 K. C0AI2O4 {x ~ 0.0) shows very little irreversibility 
and only a gentle downturn near 12 K. Below 4K, there 
is a small upturn in the susceptibility which could arise 
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TABLE I: Magnetic data of different spinel compounds. Ocw is obtained from the fit to the high temperatme inverse 
suscpetibility data as described in the text. Tjv is taken as the point where the ZFC and PC magnetization curves diverge with 
respect to temperature. / is obtained from Qcw/Tn- The inversion parameter 5 is also indicated for all samples. 
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FIG. 6: Field-cooled (FC, solid) and zero-field cooled (ZFC, 
dashed) molar DC susceptibility of the spinel compounds at 
low temperatures. Data were acquired under a magnetic field 
of 1000 Oe. 



from uncompensated spins. All other compounds with 
a; > 0.4 show the characteristic cusps of glassy systems 
associated with freezing of spins and no long range or- 
der as may be expected for crystographically disordered 
antiferromagnets. For all samples T/v was taken as the 
point of splitting between the ZFC and FC curves. 

A gradual opening of the M — H traces as x is in- 
creased is seen in Fig. [7] which accompanies the increas- 
ing concentration of Co atoms on the octahedral site. As 



the magnetic ions enter the B site the A-B interaction 
begins to dominat^^^I^, and even though both sites have 
the same number of spins, this interaction can give rise 
to uncompensated spins which could open the M — H 
loops due to the unequal number of A and B sites. The 
near neighbor antiferromagnetic interaction between Co 
on the A and B sites, which are only separated by a 
single O atom will compete strongly with the pure A- 
A interactions where the magnetic ions are separated by 
— O— B— O— linkages. This interaction thereby adds an- 
other competing exchange pathways which may prevent 
long range antiferromagnetic order between Co atoms on 
the tetrahdcral sites from being achieved and give rise 
to glassy behavior instead.'^ It should also be noted for 
high enough concentrations of Co^"*" on the B site that a 
ferromagnetic direct exchange between neighboring B-B 
atoms will begin to arise which could also give rise to the 
open M - H loops .1221 

For higher concentrations of Ga, it can be seen in Table 
[l] and Fig. [s] (a) that Qcw gradually decreases. Com- 
paring the title compounds with systems that have a well 
ordered magnetic A lattice such as C03O4 and CoRh204 
shows an almost linear dependence of Qcw with respect 
to nearest neighbor separation cIaa (and correspondingly 
next-near neighbor A- A separation). We note that the 
deviation from the linear trend for samples with x > 0.8 
can be understood by the site mixing in the samples re- 
ducing the number of magnetic cations on the A site. 
Further evidence that this the deviation from linearity is 
a result of site disorder is the fact that a separate report 
on a sample CoGa204 with half of the site mixing present 
in our samples lays on the line. 

The structural changes also affect the temperature 
where the system transitions from the paramagnetic to 
glassy state which we take to occur where the zero field 
coohng and field cooling data deviate as discussed ear- 
lier (Fig. [s] (a)). For small concentrations of Ga with 
X — 0.4, a sharp drop of the transition temperature to 
5K from 10 K for the pure C0AI2O4 which we attribute 
to the sudden increase in atomic disorder and dilution 
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FIG. 7: Magnetization as a function of magnetic field for all 
samples obtained at 2 K in fields up to 5 T. Note that each 
curve has been shifted in 3 T increments to show more clearly 
the opening of the loops with increasing x. 

of the magnetic A site lattic^^. This decrease is then 
followed by a gradual increase as more Al is replaced by 
Ga. Tristan et al. have studied the effect of replacing 
the Al^+ with non-magnetic octahedral Co'^+.^^ It is in- 
teresting to note that changes in the non-magnetic B site 
cation result in little to no change in &cw, however a 
clear increase in the ordering temperature is observed. 

Figure|8](b) plots the measured frustration parameter, 
/ = Qcw/Tn, as a function of (Iaa, the separation be- 
tween A ions in the different structures, including values 
for C03O4 and CoRh204 taken from the literature. Con- 
sidering the changes in Qcw and Tn we find with the 
exception of a; = 0.4 the frustration index decreases sys- 
tematically with increasing (Iaa separation. 

In order to better understand the behavior of the sam- 
ples, we performed classical Monte Carlo simulations of 
the J1-J2 model on the A-site spinel lattice, introduc- 
ing a certain amount of inversion randomly each time a 
simulation run was performed. We also assumed that im- 
purity spins occupying B sites interact with their nearest 
neighbors on occupied A sites via a Heisenberg exchange 
coupling Ji. 

In Fig. [9j we show the results of a scan of various pa- 
rameters. J2/J1 was set to 0.1, 0.2, and 0.3 while the in- 
version S set to 0.0, 0.05, 0.10 and 0.20. The ratio J,/Ji 
was fixed to 1.7 throughout the scan. After plotting the 
resulting inverse susceptibilities in the manner described 
in eqn. [3] we can make a few qualitative observations. 
As expected, increasing J2/J1 increases the frustration 
of the system in the sense that it lowers the temperature 
at which evidence of ordering appears. Additionally upon 
introducing even a small amount of inversion, the behav- 
ior of the susceptibility below the Curie- Weiss tempera- 
ture rapidly changes from antiferromagnetic (sharp up- 
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FIG. 8: (a) Curie- Weiss theta (Qcw) and ordering tempera- 
ture (Tjv) as a function of near-neighbor spacing between the 
A site atoms. The ordering temperature is taken as the point 
of deviation between the field cooling and zero field cooling 
data shown in figure [6] (b) Frustration index / — Qcw/Tn 
for the title compounds as a function of near-neighbor spac- 
ing between the A site atoms {cLaa increases with increasing 
Ga content, x). Also displayed is published data for C0AI2O4 
with a difi'erent degrees of inversion: 5 — 0.04. By compar- 
ison, the degree of inversion of the C0AI2O4 sample mea- 
sured in this work is close to 0.09. Values for C03O4 and 
CoRh2 04 were also obtained from reference 1331 while values 
for CoGa204 were taken from 1^21 Note that the values from 
reference |33] were not reported with lattice constants so the 
near-neighbor separation was taken from separate structure 
reports. 

ward kink) to ferrimagnetic (smooth downturn) in rough 
qualitative agreement with the experimental results. 

We also attempted to use the model to fit the indi- 
vidual susceptibility curves of the various experimental 
samples with mixed success. We first conducted simula- 
tions of the model with no inversion, in which case the 
fit was controlled by only one parameter, namely J2/ Ji- 
These simulations were all performed on a system con- 
sisting of a cube of 64 conventional unit cells (512 spins), 
though for certain values of J2/J1, simulations were also 
done using a system of 125 unit cells (1000 spins) to check 
for finite size effects, which were found to be negligible 
for the susceptibility data. 

Although all of the inverse susceptibility curves pro- 
duced from the numerics for the (5 = 0.0 case were sys- 
tematically above the experimental data, they showed 
the correct qualitative behavior in that they exhibited 
an upward turning kink, presumably at the ordering tem- 
perature Tat. By increasing the value of J2 until J2/J1 
is around 0.15 to 0.18 (and correspondigly a Ji between 
12.1 and 12.8 Kelvin), the transition temperature Tjy was 
brought into approximate agreement with the experimen- 



8 



0.8 
0.4 
0.0 
-0.4 

0.8 
0.4 
0.0 
-0.4 



a x=0.0(exp.) 
□ x=0.4(exp.) 




• 


8 = 


0.00 


* 


8 = 


0.05 




8 = 


0.10 


• 




8 = 


0.20 



(a) JJJ^ = 0.1 




0.4 0.6 
7'/|0owl 



FIG. 9: (Color online) Magnetic susceptibility results from 
Monte Carlo simulations plotted in the form of equation |3] 
compared to experimental data. From top to bottom one can 
see the change in the shape of the low temperature suscep- 
tibility curves for increasing competition between J2 and Ji. 
Within each panel the effect of atomic disorder is shown. Note 
that Ji/ Ji is held constant at 1.7 for all of these simulations. 
The dotted line represents the expected result for a perfect 
Curie- Weiss system. 



tally observed value as shown in Fig. 10 We also consid- 
ered simulations with small amounts of inversion between 
1 and 2% given that even pure C0AI2O4 exhibits a small 
amount of disorder. However, unless we ran our simula- 
tions with a value of Ji / Ji far from the range that gave 
a good fit to the doped sample with a; = 0.4 (discussed 
next), such a small amount of disorder changed the re- 
sulting susceptibilities very little. 

We have also fit the x — OA data by running simula- 
tions with an inversion of (5 = 0.36. To improve the disor- 
der averaging, the simulations were done using a 5 x 5 x 5 
supercell (1000 spin) with the data presented here being 
the average of three independent runs. As in the case 
above, all of the inverse susceptibility data produced by 
the Monte Carlo simulations was systematically larger 
than the experimental curve excluding the region below 
the downturn. Interestingly the data could be fit equally 
well by a range of J2 / Ji and Ji / Ji values, as long as both 
parameters were increased together, as can be seen from 
Fig. [To] Since it is probably reasonable to assume that 
J2/J1 decreases slightly from the value obtained for the 
X = 0.0 sample, we can take J2/J1 between 0.10 and 0.15 
and find that the ratio Ji/Ji should lay between 2.0 and 
2.3 whereas Ji is found to be between 8.5 and 9.7 K. 
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FIG. 10: (Color online) (a) Susceptibility curves obtained 
from Monte Carlo simulations where the effect of varying the 
ratio J2/J1 is tested in a system with no atomic disorder, 
(b) Susceptibility curves where the effect of varying the ratio 
J2/J1 is tested on a system with a fixed amount of atomic 
disorder around 5 — 0.36. 



It is worth noting at this point that if we judge the 
quality of the fit in a least-squares sense, the Curie- Weiss 
law describes the data better than the numerical suscep- 
tibility data given that the numerical data lies above the 
experimental data through a wide temperature range for 
which a straight line fits the experimental data almost 
perfectly. Despite this fact, it is valuable to see that the 
simulation data exhibits the same qualitative low tem- 
perature behavior as the experiments. One possible ex- 
planation for the deviation between the theory and ex- 
perimental data, aside from the assumption that only 
nearest and next nearest neighbor exchange play a role, 
is the possible over simplification of the g factor. One 
way to improve our simulations would be to allow for 
a temperature dependent g or also assigning a different 
value of g to each site. 

Finally, we attempted to fit the experimental suscepti- 
bilities for the samples with a; > 0.8 and were unable to 
obtain a good fit without changing the J2/J1 and Ji/ Ji 
by an unreasonable amount. For every set of parame- 
ters considered, the numerically calculated inverse sus- 
ceptibilities exhibited strong downturns and dropped far 
below the experimental curve after initially matching at 
high temperature. Though the reasons for this negative 
result are somewhat unclear, it seems that either the ap- 
propriate J2 / Ji and Ji / Ji ratios are very far from those 
found above or more likely the use of only two parameters 
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is an over simplification considering that atomic disor- 
der will locally modify exchange pathways and make the 
true exchange couplings site dependent. This possibil- 
ity is supported by the PDFs presented which show the 
local bond lengths of the end members is retained upon 
substitution. Such behavior may indicate that substitu- 
tion actually generates new competing J values which are 
not accounted for in our Monte-Carlo simulations rather 
than a simple modification to the existing pathways as 
may have been expected. 

V. SUMMARY 

We have attempted to understand the nature of 
the magnetic frustration in the A-site magnetic spinels 
CoAl2_xGaa;04 by substituting Ga for Al, in the hope of 
decreasing the relative magnetic coupling between near 
and next-near A atom neighbors. We have found, how- 
ever, that in the compounds presented there is a signifi- 
cant mixing of the A and B sites in the Ga-rich samples, 
and this inversion has a significant influence on the mag- 
netic coupling. The complexity of the structural changes 
which occur with substitution of Ga make isolating the 
influence of lattice expansion and site mixing a signifl- 
cant challenge. We have used Monte Carlo calculations 



to demonstrate the importance of site mixing through 
a change of shape in the simulated susceptibility curves 
which agrees closely with our experimental flndings. By 
comparing with samples that are perfectly ordered such 
as C03O4 and CoRh204, we have demonstrated that 
there appears to be a trend in which the frustration index 
/ depends weakly upon the separation between magnetic 
ions. We have also used a variety of structural charac- 
terizations and magnetic measurements to demonstrate 
how both lattice expansion and site mixing can simulta- 
neoulsy influence the frustration parameter. 
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